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Abstract. The principal aim of observational helioseismology is to determine mode parameters of the solar 
oscillations as accurately as possible. Yet estimates of the mode parameters are subject to many sources of noise, 
including inevitable gaps in data, stochastic nature of excitation processes, interferences among modes. Another 
uncertainty in frequency determination results from a modulation of the oscillation frequency in itself. It is well 
known that variation in the mean strength of the solar magnetic field modulates the frequency. We consider effects 
of the solar magnetic field variation on the spectral behavior of the power spectrum using a simple model. We show 
that the solar magnetic field variation can cause unwanted sidelobes beside the main peak in the power spectrum, 
and consequently show that the effect of the frequency modulation due to the solar magnetic field variation can 
lead a bias in frequency determinations. This effect should be considered seriously particularly when the / 7^ 
p-mode multiplets are analyzed to measure the rotational splitting since the separations in both cases could be 
comparable. In addition to a bias in frequency estimates the line width is also likely to be overestimated due to a 
line broadening, which is a measure of the life time of the mode. We therefore suggest that determination of the 
mode parameters should be derived with due care considering the magnetic field variation when a long data set 
is analyzed. The frequency modulation should be taken into account in the analysis to make a solid conclusion 
on any subtle dependence of the mode parameters, such as, that on the solar cycle. We emphasize, however, 
that a realistic model of the magnetic activity variation over the entire surface of the Sun should be involved in 
order to predict the effect of the solar magnetic field variation on the spectral profiles of the solar oscillations at 
the quantitative level. Finally, we conclude by pointing out that a new method is required to accommodate the 
stochastic force and the phase variation. 
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1. Introduction 

Helioseismology is a study of the solar oscillations sound- 
ing the internal structure of the Sun (e.g., Deubner & 
Gough 1984). Frequencies of the solar p- modes provide 
the most fundamental and important information on the 
Sun in that one may infer the solar internal structure 
and rotation from the observed frequencies (Gough 1984; 
Christensen-Dalsgaard et al. 1985). To obtain accurate fre- 
quencies, it is conventional to take a Fourier power spec- 
trum of time series data obtained from velocity/luminosity 
variation observations and to measure the frequencies of 
the peaks in the spectrum by fitting a symmetric Lorentz 
profile to the observed power spectrum, assuming that 
those frequencies are constant over the time interval en- 
compassed by the spectrum. This method also yields esti- 
mates of the line width and amplitude of the solar oscilla- 
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tions, which include clues of their damping and excitation 
mechanisms. 

The observed power spectrum of the solar oscillations 
is, however, far from an ensemble of smooth curves. Gaps 
in the time series, such as, due to the diurnal rotation of 
the Earth, produce deleterious sidelobes in the power spec- 
trum. The noisy power spectrum requires elaborate data 
reduction techniques to compensate to some extent for ef- 
fects of gaps in the data (Brown & Christensen-Dalsgaard 
1990; Lazrek & Hill 1993; Chang & Gough 1995; Fossat 
et al. 1999; Fierry Fraillon & Appourchaux 2001). Even if 
continuous data sets without gaps exist, there is still an 
uncertainty in determining mode parameters due to the 
stochastic nature of the excitation process of the solar os- 
cillations (Goldreich & Keeley 1977; Goldreich & Kumar 
1988). Woodard (1984) showed empirically that for the 
case of a harmonic oscillator excited by random noise the 
power spectrum will be distributed as the distribu- 
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tion with two degrees of freedom (for rigorous discussion 
see Gabriel 1994). That is, the standard deviation of the 
power at a certain frequency is equal to the mean power at 
that frequency. What it implies is that even if the length of 
an uninterrupted data set becomes extended the observed 
power spectrum may not converge to a smooth curve. 
Instead the observed power spectrum is likely to be more 
spiky while the resolution is improved. Alternatively, one 
may attempt to measure the mode parameters by dividing 
the whole data set into several shorter data sets and av- 
eraging obtained spectra before the fit, since the statistics 
may become more nearly Gaussian (e.g., Sorensen 1988). 

Unfortunately, however, this could not be a solution 
in analyzing the real solar power spectrum. The solar 
structure undergoes a slow change, causing variation in 
the frequency. The average of power spectra is simply an 
average of spectra at different epochs with different fre- 
quencies. Therefore, strictly speaking, the obtained fre- 
quency from the averaged power spectrum is an aver- 
age of varied frequencies, and the measured line width is 
likely to be over-estimated due to such a frequency vari- 
ation. Woodard & Noyes (1985) first reported a signifi- 
cant decrease in frequency using solar intensity data from 
ACRIM (Active Cavity Radiometer Irradiance Monitor), 
which were obtained from 1980 (near solar maximum) to 
1984 (near solar minimum). Since then the correlation be- 
tween the frequency shift and the solar cycle is firmly 
established (Libbrechtfe Woodard 1990; Woodard et al. 
1991; Bachmann & Brown 1993; Elsworth et al. 1994; 
Regulo et al. 1994; Chaplin et al. 1998; Bhatnagar, Jain, & 
Tripathy 1999; Howe, Komm, & Hill 1999; Jain, Tripathy, 
& Bhatnagar 2000; Chaphn et al. 2001). Careful inter- 
pretation of frequency dependence on the variation con- 
clusively suggests that the frequency variation is mainly 
due to the perturbation of the sound speed near the so- 
lar surface (Goldreich et al. 1991; Balmforth, Cough, & 
Merryfield 1996; Dziembowski, Coode, & Schou 2001). 

In this paper, we demonstrate that the variation of 
the solar magnetic field strength which leads to the fre- 
quency modulations of the solar acoustic modes may cause 
extra noisy peaks in the observed power spectrum in ad- 
dition to the effects of stochastic excitation and a line to 
broaden. A general introduction of the frequency modu- 
lation can be found in radio communication engineering 
(e.g., Panter 1965). In helioseismology context the solar 
magnetic field variation and g-modes were considered as 
an agent of the frequency modulation (Kennedy, Jefferies, 
& Hill 1993; Chang 1996; Lou 2001). Previous studies have 
been concentrated on effects due to the existence of g- 
modes or issues on a time series analysis. We here consider 
effects of the solar magnetic field variation on the spectral 
behavior of the power spectrum using a model and also 
show that the effect of the frequency modulation due to 
the solar magnetic field variation depends on the solar cy- 
cle, as expected in the observed correlation between the 
frequency shift and the solar cycle. The paper begins with 
a description of the frequency modulation and presents 
results of our simulations in Sect. 2. We discuss the effects 



of the solar magnetic variation on the observed spectrum 
and conclude in Sect. 3. 

2. Frequency-Modulated Solar Acoustic Modes 

As Lou (2001) adopted, one may derive the solar p- 
mode frequency shift due to the solar magnetic mod- 
ulation in terms of the variational principle formalism 
(Chandrasekhar 1964). The instantaneous angular fre- 
quency ioit) of a particular p-mode modulated by the solar 
magnetic field variation can be given as Lu(t) — luo + e/(i)> 
where ujq is the unmodulated angular frequency]^, e is the 
strength of the deviation which is assumed to be related 
to mode properties, such as, inertia, and f{t) is a function 
which is assumed to vary with the solar magnetic field 
strength. We note that e is so small that the modulus of 
loq is much greater than that of ef{t), and f{t) is a slowly 
varying function compared to the p-mode and same for all 
the solar p-modes since all the solar oscillations contribut- 
ing to the signal are experiencing the same perturbation 
agent to the acoustic cavity. We used solar magnetic field 
data|^ available at the National Solar Observatory at Kitt 
Peak to generate f{t), converting the observed magnetic 
field strength to the frequency variation with the linear 
relation obtained by Woodard et al (1991). Though the 
apparent periodicity is most likely related to presence of 
active regions that corotate with the Sun, we presume 
that this feature simulate intrinsic periodic variations of 
solar magnetic field strengths that may lead to frequency 
modulations of solar p-modes. It should be pointed out, 
however, that one should employ a more detailed model 
for the distribution of magnetic field in the solar inte- 
rior to estimate actual biases of frequency and line width 
determinations. The relation between the magnetic field 
strength and the frequency shift may not be an optimal 
function (Bhatnagar et al. 1999), but any subtle discrep- 
ancy should not be a serious problem in this demonstra- 
tion. Since Woodard et al (1991) normalized the relation 
to Z = mode of « 3 mHz, in our simulation e may vary a 
factor of two depending on ajQ we take. The strength of the 
deviation e is roughly an increasing function with the fre- 
quency. The simulated signal of the frequency-modulated 
p-mode can be written as 

huit) - A{t) cos(cjo< + e / /(r)dr -I- Vo), (1) 

Jo 

where A{t) is an amplitude which we consider to be de- 
termined such that the energy distribution follows the 
Boltzmann distribution, and ipo is a phase when t = 0. 
Without loss of generality the phase constant can arbi- 
trarily be set equal to zero and A{t) constant Aq. 

Since an arbitrary continuous function f(t) can be ex- 
pressed in terms of harmonics, according to Fourier's the- 
orem, we begin with a single-tone sinusoid as an illus- 

^ an angular frequency tx;o is 2 tt times a cyclic frequen cy uo 
^ Ftp : //argo.tuc.noao.edu/kpvt/daily/stats/mag.dat| 
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Fig. 1. Simulated power spectrum of the frequency- 
modulated mode is shown in the upper panel. The spec- 
tral line is significantly broadened and apparent sidelobes 

arc spaced at ^ 0.4 //Hz. Note that power is arbitrar- 
ily normalized and noise-free. The frequency modulation 
function f{t) is generated from the magnetic field strength 
observed in 1991 (near solar maximum) which is shown in 
the lower panel. 

tration, such that f{t) = coscUmt- Then the frequency- 
modulated function is given by 



huit) = Ao cos{u>ot + P smu>rnt), 



(2) 



where /3 is e/um so that the maximum phase deviation is 
inversely proportional to the frequency of the modulating 
function. The representation for iFuit) can be expanded 
due to Jacobi's expansion in a series of Bessel coefficients 
using the Fourier series expansion (Watson 1922). Using 
a property of Bessel functions we have 



Ipuit) = Ao ^ Jn{l3) cos{ujot + nujmt), 



(3) 



where Jn{0) is the Bessel function of the first kind with 
JniP) = {—^)^ J -n{fi) ■ Thus we have a time function 
consisting of a main unmodulated function and an infi- 
nite number of sidebands, whose relative amplitudes are 
proportional to Jn{fi) spaced at frequencies ±.nu>m- As 
[3 ^ oo, the number of sidebands inc;reases and the spec- 
tral components become more and more confined between 
Wo ± e- And the maximum height is also reduced because 
energy tends to be distributed among many small peaks. 
Note that (5 is proportional to the sensitivity of the mode 
on the solar magnetic variation. In reality, the modulating 
function is multitones in nature. The signal suffers more 
interference, causing more sidebands. Obviously the power 



Fig. 2. Similar plots with Figure 1, except that the fre- 
quency modulation function is derived from the data ob- 
tained in 1997 (near solar minimum). Sidelobes are less 
significant than those in Figure 1, yet the spectral line is 
still broadened. 

spectrum of a modulated function depends on a shape of 

the modulation function. For instance, in the case where 
the frequency modulation is a square wave function, when 
the amplitude of the square wave Aw is much larger than 
the fundamental frequency of the square wave the power 
spectrum shows two distinctive peaks in the vicinity of 
wq ± Aw. 

In Figures 1 and 2 we show the Fourier power spectra 
of simulated data generated by Equation (1). In our sim- 
ulations we presume I = mode with its cyclic frequency 
3574.7 /iHz and set e = 2 as the effect of the frequency- 
modulation should be more serious in higher frequency 
modes. The duration of the observation corresponds to 8 
months. As mentioned earlier, we take the observed mag- 
netic field strength and translate it into f{t) according 
to the linear relation given by Woodard et al (1991). We 
show two specific cases of different epochs, that is, year 
1991 (near solar maximum) and 1997 (near solar mini- 
mum) to demonstrate its dependence on the solar cycle. 
The observed magnetic field strengths we adopt have been 
shown in the lower panels separately in Figures 1 and 2. 
Besides the obvious frequency shifts, which can be defined 
such that uj[t)dt — wo where T is the interval of 

observation, unexpected sidelobes appear (Figure 1). In 
the case of year 1991, the solar magnetic modulation is 
almost a single-tone sinusoid with ~ 0.4 /iHz, which 
is a similar value of the amplitude of e/(t) in Equation 
(2). In other words, the peak phase deviation /? is ~ 1. 
This fact results in sidelobes located at ~ 0.4 /iHz apart 
from the main peak. The effect of the frequency modu- 
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Fig. 3. Similar power spectrum shown in Figure 1, except 
that all other noisy structures further than ± 0.2 /xHz from 
the main peak are set to zero. 

lation is likely more conspicuous in the case of the solar 
maximum. It should be noted, as we mentioned above, 
that these specific sidelobes at ~ 0.4 /xHz may not be 
manifest in the real solar power spectrum since the mag- 
netic variation used here with periods of ~ 30 days results 
from the apparent periodicity due to the presence of ac- 
tive regions corotating with the Sun. In Figure; 2. the only 
possibly detectable effect of the frequency modulation is a 
broadening of the main peak. In this particular case, the 
spectral line is more broadened by a factor of two than 
that in the case where there is no frequency-modulation. 
One important thing to bear in mind is, therefore, that 
according to the simulated power spectrum the spectral 
line width is likely to be overestimated in any case if the 
effect of the frequency modulation is ignored. Moreover, 
when the mode is frequency-modulated the height of the 
central peak is also likely reduced since the energy of the 
mode tends to be distributed among many small sidelobes 
as shown above. We note that a decrease in the strength of 
the modes from solar minimum to maximum is reported 
by the BiSON group, where the effect of the frequency 
modulation has been ignored in the analysis (Elsworth 
et al. 1993). The frequency modulation should be taken 
into account to make a solid conclusion on this matter. 
Otherwise, estimates of the life time and the strength of 
the solar oscillations seems likely to be biased. 

3. Discussion and Conclusion 

Magnetic activity changes the outer layers of the solar en- 
velope, modifying the resonant properties of the p-mode 
cavity and modulating the oscillation frequencies. By sim- 



ply measuring the positions of peaks in power spectra, one 
is restricted by the uncertainty in frequency determination 
due to the phase wandering induced by the changing cav- 
ity. We demonstrate that the frequency modulation by the 
solar magnetic field variation may affect estimates of the 
mode parameters, particularly during a period near the 
solar maximum. This effect should be seriously taken into 
consideration when the / ^ p-mode multiplets are ana- 
lyzed to measure the rotational splitting. As shown above 
the sidelobes due to the solar magnetic variation may be 
displaced next to the main peak at similar amounts of 
the rotationally splitting separation if the magnetic field 
varies with a roughly 30-day timescale. Although many 
attempts to correlate the solar rotation curve with the so- 
lar cycle have been made (Jimenez et al. 1994; Chaplin et 
al. 1996; Antia & Basu 2000), we suggest that conclusions 
of the dependence of the solar rotation curve on the so- 
lar cycle should be derived with due care since the mode 
estimatation could be biased due to the frequency modu- 
lation as well as the mode interference (Appourchaux et 
al. 2000). In addition to a bias in frequency estimates, the 
line width is also likely to be overestimated. In Figure 3, 
we show the power spectrum shown in Figure 1, in which 
all noisy structures further than ±0.2 /iHz from the main 
peak are set to zero, in order to emphasize the line broad- 
ening and to suppress the possible rotational effect. 

We conclude by pointing out that there was an at- 
tempt for developing a method to measure the frequency 
by fitting the temporal signal to a physical model deduced 
by the superposition of many modes, rather than using its 
power spectrum (Chang 1996; Chang & Gough 1995). This 
technique uses the fact that the unknown components to 
the wandering of amplitudes and phases of constituent os- 
cillators arising from the temporal variation of the cavity 
in which they are confined are related to each other in 
a known way. Provided that a sophisticated algorithm is 
implemented to accommodate the stochastic force and the 
surface variation causing the phase modulation, such an 
idea may yield more accurate estimates of the frequency 
and the life time of the solar short period oscillations. 
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